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Abstract 

This article! reviews some recent developments in Astroparticle Physics. Due to 
the extension of the field only part of the results and developments can be covered. 
The status of the search for Dark Matter, some recent results on Cosmic Rays and 
Gamma Ray Astronomy and the status of Neutrino Astronomy are presented. 
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1. Introduction 

Astroparticle Physics is rapidly growing into a di- 
versified field of observational and phcnomcnologi- 
cal physics, addressing questions ranging from the 
nature and distribution of Dark Matter, the mass of 
neutrinos and the origin of the Cosmic Rays, to the 
existence of antimatter. In addition the level and 
source composition of cosmologically important ex- 
tragalactic diffuse radiation levels are probed and 
tests are performed on non-Lorentz invariant inter- 
action and propagation terms predicted, e.g., in sev- 
eral ansatze for Quantum Gravity theories. 

Among the results and developments not 
covered here are the latest developments in 
gravitational wave antennas (for a recent review 
see, e^., [Q), the results on atmospheric neutrinos 
(see M), the measurements of the Cosmic Rays 
composition and spectra around the 'knee' in the 
all-particle spectrum (see, e.g., 0), or the results 
from the first flight of the AMS detector [Q. 

The structure of this paper is as follows: 
section ^ reviews the status of the Dark 
Matter search. The importance of Cosmic 
Ray measurements at low and high energies are 
discussed in section ^ followed by some of the latest 
developments in Gamma and Neutrino Astronomy 
in sections ^ andB|, respectively. A brief outlook is 
given in section 

2. The Dark Problems 

Some recent reviews on cosmology, e.g., [^, 
conclude that for the first time now we have a 
complete observational accounting of matter and 
energy in the Universe consistent with ACDM 
(i.e.. Cold Dark Matter plus non-zero A term) 

t Invited review at EPS-HEP '99, July 1999, Tampere, 
Finland. 



inflationary cosmology. This, however, is not 
equivalent to a complete understanding of these 
ingredients, and according to one summary shown 
in figure H the three Dark Problems that we are 
facing today are: Where are the dark baryons? 
What is the non-baryonic Dark Matter? What is 
the Dark Energy recently detected in analyses of 
distant supernovae In the following the status 

of the questions concerning the amount and nature 
of Dark Matter (DM) will be addressed. 

The evidence for DM is manifold and can 
roughly be divided into evidence on galactic scales 
(i.e., 0(10) kpc) and on galaxy cluster scales (i.e., 
0(10) Mpc). Measurements of rotation curves of 
galaxies js), especially of spiral galaxies including 
the Milky Way, still give the most convincing 
evidence for the presence of a galactic Dark Matter 
(DM) mass component which is about a factor 10 
larger than the observed luminous matter. Most 
of this mass has to be distributed in an extended 
halo in an approximately spherical distribution with 
density falling as 1 /r^ . Inspection of figure |l| shows 
that this factor can in principle be accommodated 
by baryons produced in Big Bang Nucleosynthesis 
(BBN). On scales of the extension of rich galaxy 
clusters the second DM problem becomes visible. 
The main mass fraction of rich clusters is in form 
of hot intracluster gas, the mass of which can be 
inferred by measuring the X-ray flux or by mapping 
the Sunyaev-Zel'dovich decrement caused by the 
scattering of cosmic microwave background (CMB) 
photons on the hot electrons. The total cluster mass 
is then derived by one or more of three methods: 
from motion of the galaxies utilizing the virial 
theorem, by assuming hydrostatic equilibrium for 
the hot gas, or, in some cases, by gravitational 
lensing. From these data the mean matter density 
in the Universe, f^M, can be inferred. The obtained 
value of r^M = 0.4 ± 0.1 is consistent with values 
derived from other observables like bulk matter 
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Figure 1. Current status of the observational and 
theoretical census of matter and energy in the Universe. 
From iBl. 



flow or the peculiar motion of our Local Group of 
galaxies (e.g., |^). Assuming BBN to be correct 
the main matter component in the Universe thus 
has to be non-baryonic. Note that the validity of 
BBN is well tested by the measured ''He, D, '^He, 
and ^Li abundances (e.g., |^). In addition to the 
data there are solid theoretical arguments for the 
existence of a large non-baryonic DM mass fraction 
in the Universe. E.g., without large amounts of 
non-baryonic DM there is no working theory of 
galaxy formation consistent with the small level of 
fluctuations observed in the CMB. 

So far DM has only be discovered by its 
gravitational effects and the search for the 
manifestation of DM in other interactions mainly is 
limited to searches for galactic DM. These searches 
can either focus on the baryonic or a possible 
non-baryonic component. Note that a debate on 
whether any non-baryonic component is present at 
galactic scales at all has been going on since the first 
solid evidence for galactic DM was obtained. 



Candidates for the baryonic component are low 
mass stars, stellar remnants, Massive Astrophysical 
Compact Halo Objects (MACHOs) such as brown 
dwarfs and Jupiters with masses below ^ 0.09 
M0, cold molecular clouds, or very massive black 
holes. Of these candidates all but the cold molecular 
clouds are in principle detectable by gravitational 
microlensing. 

Gravitational microlensing (ML) as a means to 
search for MACHOs was proposed in a seminal 
paper by Pacyhski in 1986. Only 13 years 
after this proposal ML is a thriving field and can 
be considered a new field of galactic astronomy 
(for a review see, e.g., ||ll|). Several different 
target directions have been monitored in the search 
for the transient ML events characterized by the 
achromatic, symmetric and unique amplification of 
a source star. The pioneering work was performed 
by the EROS and MACHO collaborations who 
in the last years monitored millions of stars in 
the small (SMC) and large (LMC) Magellanic 
Clouds (see figure ^) and announced the first 
candidate events in 1993 jlj, |3). The total 
ML events statistics towards the SMC and 
LMC currently is about 20 events (from EROS, 
MACHO, and OGLE). Other target directions 
have been the galactic bulge (e.g., OGLE, OGLE- 
2, DUO, MOA, and MACHO), galactic spiral 
arms (e.g., EROS), or towards the Andromeda 
galaxy (AGAPE). In addition new collaborations 
(e.g., MACHO/GMAN, MPS, or PLANET) aim at 
accurate photometry of ML events in order to search 
for planets in the source star systems. The ML 
event statistics towards galactic targets is currently 
more than 400 and is utilized to study details of 
galactic structure unobservable by other methods 
(see, e.g., Q). 

For the interpretation of ML events in terms 
of DM one has to take into consideration that the 
velocity, distance and MACHO mass degeneracy in 
the basics ML equations (see Eq. (1) and (2)) does 
not allow to determine the location of the lenses 
itself based on the magnification and duration of 
the event alone. Only for special events like binary 
ML events discovered towards the Magellanic clouds 
and long duration events showing parallax effects 
due to the movement of the earth around the 
sun is it possible to constrain the location of the 
lenses. E.g., for the binary-lens caustic crossing 
event MACHO-98-SMC-1 which was monitored by 
5 ML collaborations after the MACHO alert, the 
combined analysis yields strong evidence for the lens 
system to be located in the SMC itself Q . 
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Figure 2. Schematic view of lensing populations 
discussed for the Galaxy: the standard spherical 
halo (solid lines), a heavy spheroid (dashed lines), a 
maximum thick disk (solid lines) and a dark thin disk 
(dotted lines) . For each population the density contours 
p = 10~^ and p = 10"^ Mo/pc^ are shown. The 
locations of the sun and the small (SMC) and large 
(LMC) Magellanic Clouds are indicated. From 



In general, however, the basic ML equations 
only relate the velocity, distance and mass of the 
lenses to give the lensing probability: 



T oc p{x)x{l — x)dx 



(1) 



with the mass density distribution, p(x), of 
microlensing matter needed as an input according 
to the adopted isotropic halo model (see figure 

and X = Dlens I D source with Di^ns and D source 

the distances of the lens and the source star, 
respectively. The duration of the lensing event is 
given by: 



At (X 



lMx{l 



(2) 



with M the mass and v± the transverse velocity 
of the lens. The lensing probability, r, gives the 
probability for finding a source star within the 
Einstein radius, Re, of some MACHO and which 
thus is magnified by a factor A — 1.34. The 
smallness of r (O (10~^)) makes it necessary to 
simultaneously monitor millions of stars. 

The combined limits and results on the galactic 
halo mass fraction in form of MACHOs from the 
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Figure 3. Exclusion diagram at 95% CL assuming the 
standard spherical halo. The dashed line is the limit 
from EROS 1 towards the LMC and SMC ||l7l, the solid 
line is the Umit from EROS 2 towards the SMC . The 
95% CL allowed region from the MACHO collaboration 
is shown as the shaded area |l8[ with the preferred value 
indicated by the cross. The tnm dashed line corresponds 
to the limit obtained under the assumption that no halo 
event has been observed. From [hdl . 



observations of the LMC and SMC by the EROS 
and MACHO collaborations are shown in figure |^ 
p6[ . From the lack of short duration events it can 
be concluded that MACHOs in the mass range 10~^ 
Mq S TO < 10-3 Mq make up less than 25% of 
the halo mass /or most halo models. On the other 
hand, the MACHO collaboration derived a most 
probable MACHO mass range of TO = O.dtoi Mq 
from the observation of the average long duration 
of 8 candidate ML events in the first 2.1 years 
of LMC data jl^. In addition this event rate is 
compatible with about 50% of the mass of the halo 
in the form of MACHOs of this mass. In figure || 
the 95% CL allowed halo mass fraction region as a 
function of tomacho derived from these data is also 
shown. Recently the EROS collaboration published 
its results on the first 2 years of observation towards 
the SMC jlra. Only one ML candidate event was 
observed. This low event rate can be transformed 
into a limit on the fraction of the halo made out 
of heavy MACHOs assuming a standard spherical 
halo model. The 95% CL exclusion limit is shown 
as the heavy solid line in figure ^ and now excludes 
a halo completely made out of heavy (i.e., ^ 0.5 
Mq) MACHOs. Note that these limits are lower 
if some of the observed ML events are due to self- 
lensing. Self-lensing means that the lenses are not in 
the halo but in tidally elongated Magellanic clouds 
themselves or in the (warped) disk of our own 
galaxy. In addition the lensing events might be 
due to faint stars of an yet undiscovered intervening 
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dwarf galaxy. If the ML events should be caused by 
one of the above effects MACHOs will be irrelevant 
for DM whatever their astrophysical nature. In 
order to study this further and get a better handle 
on the location of the lenses through the comparison 
of mass moments extracted from the ML data with 
halo models (e.g., a much larger event statistics 
towards the LMC and SMC and other line-of-sights 
are necessary. 

Independent constraints on the DM contribu- 
tion of baryonic compact objects come from astro- 
physics and cosmology which now rule out most 
of the compact candidates for baryonic DM in the 
halo. Without discussing each limit in detail the 
situation of protostellar objects, stellar objects and 
stellar remnants as candidate objects for making up 
the halo mass can be summarized as follows: 

• Jupiters (10~'^ < M < 10"^ Mq) are ruled out 
by ML data, i.e., halo mass fraction < 25%. 

• faint hydrogen burning stars (M > 0.09 Mq) 
and young brown dwarfs (M < 0.09 Mq) are 
ruled out on the basis of HST and ISO data 
and contribute a halo mass fraction of ~ 0.01 

• old brown dwarfs (M < 0.09 M©) are ruled out 
by HST and US Naval Observatory parallax 
data and star formation theory and contribute 
a fraction 0.03 p^]. 

• white dwarfs (wdj^ncutron stars, and black 
hole stellar remnants are ruled out by 
observations strongly constraining the number 
of possible progenitor stars. The constraining 
data range from chemical abundances of C and 
N measured in Lya systems , global D and 
'''He abundances to the stringent limit on 
the cosmic infrared background flux inferred 
from the observation of multi-TeV 7-rays from 
the blazar Mrk 501 at a redshift of z = 0.034 

. All of these signal levels would have been 
raised by progenitor stars, yielding: 

- n^dh < 2xlO"M||l 

- n^dh < 0.003 (Hi 

- ^wdh < (1-3) xlO-3 

with the Hubble parameter Hq = • 100 km 
s-i Mpc-^ 

In conclusion, the ML data in combination with 
current astrophysics seem to rule out a completely 
baryonic halo, if it is made out of isotropically 
distributed compact massive objects. The nature 
of the observed events towards the LMC and SMC 
remains unknown. One possible explanation is a 
halo made out of primordial black holes with a mass 
of about one solar mass which would explain the ML 
events and act as cold DM otherwise. Other exotic 



candidates like mirror matter stars which would also 
circumvent the astrophysical limits are discussed in 
the literature (e.g., P^). An astrophysics scenario 
which would still allow for a purely baryonic halo is 
a halo made out of cold molecular clouds (see, e.g., 
psi ) without explaining the ML events. However, 
most of the data and theoretical models currently 
point to a substantial fraction of the halo being 
made up by non-baryonic DM. 

2.2. Search for Non- Baryonic Dark Matter 

Experiments have to answer the question of 
the nature of the non-baryonic DM (NBDM). 
The search for observable signals resulting from 
interactions other than gravity in general require 
detectors especially suited to the DM candidate in 
question. 

At the moment there is a surplus of matter 
(types) in cosmology and at the same time 
important stable and massive particles predicted 
to exist either within the Standard Model or 
in its plausible extensions have not yet been 
discovered. The most promising approach to 
searches for NBDM therefore lies in performing 
dedicated searches for selected Particle Physics 
candidates. These have to fulfill the requirement 
to have decoupled from the primeval plasma when 
their interaction rate became smaller than the 
cosmic expansion rate and which therefore would be 
floating around the Universe today as relic particles. 
According to their energies at the time of decoupling 
one differentiates between relativistic (for masses 
less than 1 keV), i.e., hot, and non-relativistic, 
i.e., cold, relics. The best candidates which could 
contribute significantly to the mass-energy density 
in the Universe (i.e., SI ~ 1) are the axion, the 
neutrino, the lightest supersymmetric particle if R- 
parity is conserved, or magnetic monopoles. Note 
that monopoles do not work well in numerical 
simulations of structure formation and are thus 
regarded to be the least probable solution to the 
NBDM problem by most cosmologists and we shall 
in the following not discuss experimental searches 
for monopoles. 

Of the remaining candidates the axion is a very 
good candidate but hard to detect in laboratory 
experiments within the current mass constraints. 
It is the Goldstone boson of the U(l) Peccei- 
Quinn symmetry [E8| as the elegant solution to 
the strong CP problem of QCD. At the same time 
it would be a good cold DM (CDM) candidate 
as, within the current laboratory and cosmological 
limits, relic axions could make up a significant 
fraction of the critical density. The sensitivity 
of experiments searching for light axions via the 
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PrimakofF conversion of axions into microwave 
photons in strong magnetic fields are now entering 
into the cosmological interesting range . See js^l 
for a detailed review of the recent development 
and the current limits which mainly stem from 
astrophysics and cosmology 

The next candidate, a light massive neutrino, 
would be a hot DM (HDM) particle and its mass 
would be directly proportional to its contribution to 
il. More precisely, the sum of the mass eigenstates 
which contribute to the active flavour neutrinos 
is directly proportional to fi^. In order for the 
structure formation theories to work neutrinos (or 
other HDM) cannot be responsible for the whole 
non-baryonic DM fraction but have to give a 
contribution to ft of ~0.15. This can be translated 
into J2 eV (for Hq = 65 km • s'^ • 

Mpc-i and i = 1, 3) ||. This number, 

derived from a cosmological measurement, i.e., 
the fluctuation spectrum of the cosmic microwave 
background, is compatible with (and competitive 
to) the direct determination of the mass of the 
electron antineutrino in the measurement of the 
endpoint of the (3 spectrum in tritium decay. 
These experiments now start to consistently yield 
a positive antineutrino mass in the fit, with the 
current value from the Mainz experiment reported 
at this conference of mp < 2.8 eV (95% CL) 
In spite of the large relic neutrino density (54 
cm~'^ per light flavour) the determination of their 
exact contribution to 17 will most probably come 
from this type of direct laboratory experiments, 2/3 
decay experiments (e.g., |3|]), or the forthcoming 
long-baseline oscillation experiments (e.g., [||). 
However, as the 2/3 decay experiments only give 
a mass convolution involving unkown phases and 
oscillation experiments only give mass differences, 
the improvement of the cosmological limits will 
continue to be important. 

The fundamental questions that are left open 
by the Standard Hot Big Bang Cosmology point 
towards a grander theory. The best candidate 
to date is Inflation plus cold Dark Matter 
leading to a flat and slowly moving Universe. 
In addition Particle Physics models beyond the 
Standard Model offer a variety of CDM candidate 
particles. As Superstring theories predict a 
Supersymmetry sector at the TeV scale the 
favourite CDM candidate in extensions of the 
Standard Model is the Lightest Supersymmetric 
Particle (LSP), presumably the neutralino, in 
minimal supersymmetric models with conserved 
R-parity. If these particles have been produced 
in the early Universe and have decoupled when 
they became non-relativistic their density would 
be inversely proportional to their interaction rate. 



The requirement that the LSP contributes to Q. 
of 0(1) makes it weakly interacting with a mass 
between -10 GeV and ~500 GeV [||, hence 
Weakly Interacting Massive Particle (WIMP). 

Assuming that WIMPs contribute a substantial 
fraction to the galactic halo mass leads to a density 
at the location of the earth in the range of 0.3 
- 0.4 GeV cm~'^ in the isothermal standard halo 
picture (depending slightly on the chosen halo 
parameters). WIMPs can be detected directly or 
indirectly, where their elastic scattering on matter 
constitutes the direct channel and the detection 
of annihilation products expected from pairwise 
annihilation processes define the indirect channels. 

The direct searches are based on the detection 
of nuclear recoils of O(keV) as signals of the elastic 
scattering process within the sensitive detector 
volume. The expected signal distribution is roughly 
exponential with the mean depending on the mass 
of the WIMP. The expected signals rates are of the 
order of 1 event per day and kg of detector mass. In 
the scattering only a small fraction of the delivered 
energy, however, goes into the ionization channel. 
Most of the energy is released as heat and cryogenic 
devices are therefore intensely studied for their 
suitability. The main background is due to electron 
recoils from Compton scattering of background 7- 
rays. Due to the motion of the earth through the 
WIMP halo an annual modulation of the direct 
scattering signal with an amplitude of the order of 
a few % is expected. 

Detectors with a low inherent radioactive 
background like Germanium, Nal, or cryogenic 
detectors are mostly used to search for WIMP 
scattering signals (e.g., for a review). The 
DAMA collaboration presented a hint for an 
annual modulation signal in 1997 ||35|| and recently 
announced ~3cr evidence for annual modulation 
after quadrupling the analysed statistics p6| . 
The residual detector count rate versus time 
is shown in figure g The 90% CL allowed 
region of scalar WIMP-nucleon scattering cross 
section versus WIMP mass only barely enters 
into the region expected in constrained minimal 
supersymmetric models (CMSSM) (see figure 0). 
The constrained or 'phenomenological' MSSM is the 
simplest supersymmetric extension of the Standard 
Model with the number of free parameters of the 
most general R-parity conserving model of 124 
reduced to 7 real parameters. Two other Nal 
experiments, the UK Dark Matter Collaboration 
(UKDMC) [|0| which is taking data since 1996 
in the UK Boulby Mine and a French experiment 
at Modane |Q studying Nal(Tl) crystals for their 
suitability as DM detectors, use the PMT pulse 
shapes to identify and distinguish against a, /3, 
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Figure 4. Residual count rate versus time elapsed since 
January 1 of the first year of DAMA data taking. The 
expected modulation is a cosine function with minimum 
at December 2 and the maximum at June 2. From [B7|. 
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and 7 backgrounds. Both experiments find a 
not yet understood population of short pulses not 
consistent with background expectation according 
to elaborate background studies (e.g., These 
data are also not fully compatible to the WIMP 
signal expectation and are thus not interpreted in 
terms of a WIMP signal. Due to the poor intrinsic 
background rejection capability of Nal detectors 
rather some hidden systematics is suspected. In 
order to test the DAMA evidence in future searches 
the active identification of the background will be 
of foremost importance (e.g., with liquid Xenon or 
cryogenic devices). Note that, e.g., the current 
UKDMC sensitivity is comparable to DAMA. Some 
of the existing limits on the scalar WIMP-nucleon 
cross section versus WIMP mass together with 
some expected sensitivities for some of the planned 
experiments are shown in figure ^. 

Indirect searches for WIMP DM are based on 
the detection of the products of pairwise WIMP 
annihilation in the galactic halo or in the centre of 
earth or sun. This process will occur if the WIMP is 
the supersymmetric neutralino, X: a Majorana 
fermion. One of its characteristics is the production 
of equal amounts of matter and antimatter and 
all stable annihilation products, e.g., p, p, e~, e+, 
D, u, may serve as signature. In addition 7s can 
be produced in loop-induced annihilation reactions 
Q,^. Whereas protons and electrons are too 
abundant in the normal Cosmic Rays, antimatter, 
which does not exist in sizable quantities in the 
observable Universe, may herald the annihilation 
of neutralinos in the galactic halo. The expected 
annihilation rate depends on the WIMP density 
distribution in the halo. According to the halo 
models shown in figure ^ and according to results of 
N-body simulations of the formation of dark halos 



Figure 5. WIMP-nucleon cross section limits for 
scalar interactions as a function of WIMP mass. The 
shaded area indicates the DAMA 90% CL evidence 
contour. The solid lines indicate excluded regions. The 
broken lines indicate the expected sensitivity for the 
ongoing CDMS experiment at different site s ||39| ] and 
the proposed GENIUS experiment [ [43| . See [P8|[ for full 
references. From |j8| . 

a considerable enhancement of the DM density near 
the Galactic centre (GC) is expected. For the x X 
— > 7 Xq channel this will lead to a much enhanced 
signal from the direction of the GC where Xq is 
a neutral particle of mass mx and the 7's will be 
nearly monochromatic with an energy 



Due to the non-relativistic velocities of the WIMPs 
this channel thus is characterized by a huge signal 
peak-to-width ratio with no known astrophysical 
background source. The proposal to use ground- 
based Imaging Air Cherenkov Telescopes (lACTs) 
to search for this signal dates back to 1992 [E6| 
when it was realized that the only average energy 
resolution of the order of 10-20% is offset by 
the huge collection areas of the order of 10^-10^ 
m^. According to recent calculations |^ the 
upcoming generation of lACTs (see section 4) will 
provide enough sensitivity to probe a significant 
part of the parameter space of the CMSSM. In 
figure ^ an estimate of the flux sensitivity of the 
MAGIC Telescope currently under construction on 
the Canary Island La Palma |48| is shown in 
comparison to the flux predictions for the channel 
XX ~^ 77 for a large number of CMSSM models pTf . 
In addition figure ^ shows the expected coverage 
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Figure 6. Flux sensitivity (5 ct in 50 h observation 
time) of the upcoming MAGIC Telescope (full line) for 
a 7-ray source at the position of the GC compared 
to flux calculations within GMSSMs. Each point 
represents the prediction for one possible CMSSM 
model representation. All points satisfy 0.025 < Sl^/i^ 
< 1. From 



for the CMSSM parameter space covering the 
DAMA evidence region illustrating the somewhat 
orthogonal sensitivity to the supersymmetric 
parameter space in terms of direct and annihilation 
cross section coverage for the same mass range. 
A final remark on the DM problem: recent de- 
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Figure 7. As for figure ^ but only the flux predictions 
for those GMSSM models are shown which predict a 
scalar elastic scattering cross section compatible with 
the DAMA evidence contour. From U^fi. 



velopments in experimental and phenomenological 
Particle Physics point towards the fact that the ex- 
istence of the observed baryon asymmetry itself is a 
signature of physics beyond the Standard Model. 
Csikor et al. showed in a 4-dim. lattice calcula- 
tion jsi],^ that there is no first order electroweak 
phase transition in the Standard Model if muiggs ~ 
75 GeV. The critical point in the phase diagram of 
the electroweak phase transition is thus below the 
current lower Higgs mass limit from LEP II p3| . 
Therefore one of Sacharovs conditions for produc- 
ing the observed baryon asymmetry (baryon num- 
ber violation, CP violation, and out-of-equilibrium 
dynamics as in a first order phase transition) is 
not fulfilled for the electroweak symmetry break- 
ing phase. Another mechanism which may be re- 
sponsible for the baryon asymmetry and which has 
received growing attention lately due to the exper- 
imental evidence for a non-zero neutrino mass, is 
the leptogenesis mechanism [Q. In this model 
the cosmological baryon asymmetry is generated 
from a primordial lepton asymmetry which was pro- 
duced by the out-of-equilibrium decay of heavy Ma- 
jorana neutrinos. In a leptogenesis scenario Bolz, 
Pliimacher and Buchmiiller p5| showed that with- 
out encountering the 'gravitino problem' the non- 
baryonic DM might be present in the form of gra- 
vitinos (as the lightest supersymmetric particle) 
with masses in the range between 10 GeV and 100 
GeV. This would mean that the nonbaryonic DM 
would be present as GIMPs, i.e., Gravitationally 
Interacting Massive Particles, in which case all lab- 
oratory searches would be negative and only astron- 
omy, e.g., through weak gravitational lensing and 
ML, could provide more data on the DM. 

3. Cosmic Rays and Antimatter 

In the last years advancements in detector 
technology have led to an intensified study of many 
aspects of Cosmic Ray physics. With the upcoming 
Cosmic Ray (CR) experiments like, e.g., AUGER, 
AMS, and PAMELA, and the upcoming 7-ray 
detectors, e.g., CANGAROO II, GLAST, HESS, 
MAGIC, and VERITAS, it is safe to predict that 
this trend will continue (see also section 4). As the 
'beam' of CRs is either background or signal the 
understanding of the CR spectrum and composition 
is of utmost importance for many current and 
upcoming astroparticle physics experiments. 

Here we only report on some selected results and 
a few of the upcoming experiments focusing on the 
low energy and highest energy part of the measured 
CR energy spectrum. In the low energy domain 
the interest in CR measurements, beyond the 
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Figure 8. BESS 1995 and 1997 (solar minimum) 
antiproton fluxes measured at the top of the 
atmosphere. Also shown are some of the previously 
existing data. The error bars correspond to the 
quadratic sums of the statistical and systematic errors. 
The curves represent some recent calculations of the 
expected spectra from CR interactions with interstellar 
matter for the solar minimum period. From mM. 



understanding of the sources and propagation of the 
CRs themselves, is focussed on the derivation of the 
neutrino oscillation parameters, the measurement 
of the antiproton and positron fluxes in the search 
for neutralinos and primordial black holes, and the 
search for antimatter. 

A large fraction of past and present data 
on the absolute CR flux come from balloon- 
borne experiments. A representative detector 
flown in this type of experiment is the Japanese 
BESS spectrometer. The main parameters and 
components of this detector are the 1 Tesla 
magnetic field produced by a thin (4 g/cm^) 
superconducting coil filling a tracking volume 
equipped with drift chambers providing up to 28 
hits per track with an acceptance of 0.3 m^ sr. In 
addition two hodoscopes provide dE/dx and time- 
of-flight measurements. In the 1997 flight data 
at the top-of-the-atmosphere (mean residual air 
mass 5.3 g/cm^) were collected for 57,000 s. The 
combination of the measured p fluxes as a function 
of kinetic energy from the 1995 and 1997 campaigns 
are shown in figure @ pa] . 



Antiprotons should be produced as secondaries 
in interactions of the galactic CR protons with the 
interstellar medium. Their kinetic energy spectrum 
is expected to show a characteristic peak around 2 
GeV, with sharp decreases in flux towards smaller 
and larger energies. Other possible sources of p are 
the annihilation of neutralinos at the GC or the 
evaporation of primordial black holes (PBH). The 
latter possibility received growing attention in the 
recent past following the detection of antiprotons 
with kinetic energies below 0.5 GeV by BESS in 
1993 ||5^ . Today black holes are only expected to be 
formed in stellar collapses of stars with several solar 
masses. The uniformity of space-time precludes 
collapses of matter to black holes with masses 
below this limit. In the early Universe, however, 
virulent conditions may have led to the formation 
of PBH with arbitrarily small masses |Q , e^. , by 
the collapse of large density perturbations 59, 60|. 
Data constraining PBH abundances will thus yield 
constraints on the density fluctuation spectrum 
in the early Universe, an important ingredient 
of structure formation theories. Quantum effects 
lead to the evaporation of these PBH by particle 
emission | |6l| . The emission spectrum is similar to 
a black body with finite size and a temperature 
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where Mpbh is the mass of the PBH in grams. 
For TpBH above the QCD-scale of - 300 MeV 
the evaporation process will result in relativistic 
quarks and gluons which may produce antiprotons 
during hadronisation. The expectation for the 
flux of antiprotons from this process is a spectrum 
increasing towards lower kinetic energies down to 
^^0.2 GeV |62|. Similarily the spectrum of p 
expected from the annihilation of neutralinos at 
the GC is characterized by a significant flux below 
1 GeV [Hj. The uncertainty of the expected 
secondary flux due to the uncertainty of CR 
propagation in the Galaxy, however, is still of the 
same order of magnitude as the 'signal' expected 
from neutralino annihilation. 

This propagation uncertainty is illustrated in 
figure H The data points show that BESS now 
has convincingly measured the predicted peak in 
the kinetic energy distribution at ~2 GeV. The 
location and height of this peak can be calculated 
by using the measured p and He CR spectra 
together with the cross sections for p production 
measured at accelerators. The different propagation 
model calculations shown in figure ^ differ mainly 
in propagation parameters. Although the data 
below 2 GeV now are much better from the 



statistical point-of-view, the uncertainty in the 
propagation of secondary ps still precludes any 
firm conclusion on additional 'direct' sources of 
ps. Future measurements at times with higher 
solar activity will be performed, e.g., by BESS, in 
order to study the secondary component in more 
detail. At times of high solar activity the 'direct' p 
component should be virtually absent. In addition 
extending the measurements to larger energies (i.e., 
Ekin 3> 2 GeV) will also yield constraints on the 
propagation models. Here the current experiments 
did not have enough sensitivity and exposure to 
limit the propagation parameters. Upcoming CR 
experiments like, e.g., the PAMELA experiment, 
will here yield much improved data. PAMELA 
is a satellite-borne magnet spectrometer currently 
being built by the WiZard collaboration together 
with Russia. It is planned to be launched from 
Baikonur in 2002. The scientific objectives of the 
PAMELA mission are to measure the spectra of 
p, e+, and nuclei in a wide range of energies, to 
search for primordial antimatter and to study the 
CR fluxes over half a solar cycle. PAMELA will be 
able to measure particles with magnetic rigidities 
(momentum/charge) up to a few hundred GV/c. 

Contrary to the antiproton signal the positron 
signal from neutralino annihilation at the GC is 
expected to show up at large kinetic energies, 
i.e., well above the geomagnetic cutoff. The 
expectations for improved data in this channel from 
the upcoming experiments are also high. Similar 
to the antiproton channel the current data base on 
CR positrons does not yet allow any conclusion on 
additional sources besides CR interactions in the 
interstellar medium. 

Another important aspect of CR physics is 
the search for primordial antimatter. This search 
mostly focuses on antihelium nuclei and currently 
the upper limit on the CR anti-helium/helium ratio 
provides the best evidence for the Galaxy and the 
nearby Universe being made up solely of matter. 
The combination of the BESS collaboration data 
taken during 5 campaigns between 1993 and 1998 
have led to an improvement of the upper limit on 
the ratio of the abundances of antihelium to helium 
with magnetic rigidities between 1 GV/c and 16 
GV/c by about a factor of 90 [|5|. The upper 
limit of 10~^ quoted by BESS is compatible with 
the limit presented by AMS Q and any further 
substantial improvement will need long duration 
measurements, e.g., by PAMELA or AMS. 

The absolute fluxes of proton, helium and 
atmospheric muons are important for the derivation 
of the neutrino oscillation parameters since the 
atmospheric neutrino flux is proportional to the 
normalization of the dominating CR proton and 
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Figure 9. The absolute differential CR proton 
spectrum as measured by the BESS collaboration in 
1998. Also shown are some previous measurements, 
and, as the full line, the spectrum assumed in the Honda 
et al. pm calculation of the atmospheric neutrino 
fluxes. From H. 



helium fluxes. In addition there other applications 
where the absolute measurement of the proton 
spectrum is important, e.g., in predictions of the 
secondary antiproton and positron fluxes needed to 
constrain exotic production channels and the above- 
mentioned propagation models. 

Figure shows the most recent measurement 
of the CR differential proton spectrum in the 
kinetic energy range between 1 GeV and '^100 GeV 
by BESS to geth er with data taken by previous 
experiments [|67| . Shown as the full line is the 
proton spectrum assumed in the Honda et al. 
Monte Carlo calculation used in the analysis 
of the Su per Kamiokande data on atmospheric 
neutrinos [p8|. The different data sets show 
discrepancies in the measured absolute fluxes as 
large as a factor of 2 at '^50 GeV. The absolute 
differential proton flux was also measured with the 
MASS balloon detector in 1991 but published only 
recently [|69| . In addition the differential helium flux 
at the top of the atmosphere and the differential 
muon fluxes in the atmosphere as a function of 



10 



> 

^10 = 





25.5 



10 10 10 10 10 10 

Atmospheric depth (g/cm^) Atmospheric depth (g/cm^) 



Figure 10. Differential (a) /i~ and (b) fi'^ fluxes as a 
function of atmospheric depth. From top to bottom the 
momentum ranges (in GeV/c) are: 0.3-0.53 (scaled by 
lO'^), 0.53-0.75 (10*), 0.75-0.97 (10^), 0.97-1.23 (10^), 
1.23-1.55 (10), 1.55-2 (1), 2-3.2 (1), 3.2-8 (1), and 8- 
40 (1). The /i+ data are shown up to 2 GeV/c. The 
solid lines represent the Bartol Monte Carlo calculation 
results [iTi. From frT 



atmospheric depth were measured. The proton 
and hcUum fluxes above 10 GeV were found to be 
compatible with the LEAP '87 data |Q as shown 
for the proton case in figure ^ corroborating a lower 
normalization of the absolute proton and helium 
fluxes as compared to previous compilations of data, 
e-g-, iPl- 

Tne CAPRICE94 balloon experiment measured 
the atmospheric muon flux as a function of 
atmospheric depth and compared the data 

to calculations using the Bartol Monte Carlo |7^ 
(figure nG). The disagreement between data and 
Monte Carlo prediction increases with atmospheric 
depth from about 1.1 ± 0.1 (stat.) ± 0.1 (syst.) for 
the ratio of simulation result to the measured /u" 
flux at the top of the atmosphere up to 1.8 ± 0.1 ± 
0.1 at the maximum of the muon flux at a depth of 
around 200 g/cm^. As the atmospheric muon flux is 
closely related to the flux of atmospheric neutrinos 
10 these data will have to be used to improve the 
Monte Carlo generators. Gaisser investigated the 
influence of a change in slope and normalization 
of the input proton spectrum on the atmospheric 
neutrino result |^ . He concluded that for a change 
in the flux normalization above 10 GeV from the 
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Figure 11. The top end of the CR energy spectrum 
as measured by AGASA, Fly's Eye, Havarah Park, and 
Yakutsk. The consistency of the spectrum is illustrated 
by the fact that only a shift in the absolute energy scale 
of ~15% was necessary in order to align the spectra. 
Note that the measured flux values were multiplied by 
E"^ to enhance the 'ankle' feature at around 10^^ eV. 
From ifTi. 



flux assumed by Honda et al. [|66[ in the direction of 
the new BESS result the observed excess of electron 
neutrinos would be increased by a similar amount 
and the deficit of muon neutrinos correspondingly 
reduced. Taking these new data into consideration 
in the fit of the SuperKamiokande data will very 
likely yield a shift in the fit parameters. Such a 
shift could be of great importance for the upcoming 
long-baseline neutrino experiments and illustrates 
the strengthening of the interconnection of Particle 
Physics and Astroparticle Physics in the search for 
New Physics. 

At the uppermost end of the CR energy 
spectrum in the past 40 years a number of 
experiments have been collecting events with 
energies larger than lO'^^ eV. The galactic magnetic 
field (B '-^ 3 /uG) cannot contain CR protons 
with momenta larger than a few times 10"'^* eV/c. 
CRs with momenta above this value are thus very 
likely of extragalactic origin. Figure |ll| shows a 
compilation of the all-particle flux above 10^^ eV as 
measured by the four experiments AGASA, Fly's 
Eye, Havarah Park, and Yakutsk until 1995 [[Ta . 
The highest recorded CR energy to date is 3 x 10^" 
eV [0 . Of the four experiments which contributed 
to figure |ll| the AGASA (Akeno) group has reported 
new results this year. The Akeno 20 km^ air 
shower array was operated in Japan from 1984 to 
1990 and became part of the still operational '^100 
km^ Akeno Giant Air Shower Array (AGASA) in 
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1990. The detection of an anisotropy in the arrival 
direction of CRs with energies around 10^^ eV (0.8 
- 2 X 10^8 eV) was reported by AGASA in 1998 
Q and updated recently The amplitude 

of the anisotropy in the first harmonic analysis 
was found to be 4% and could be identified in a 
2-dimensional map as stemming from 4.5 cr and 
3.9 a excesses from the direction of the Galactic 
Centre and the Cygnus region, respectively. This 
is a spectacular result as it is clear evidence for 
the existence of galactic CRs up to this very high 
energy. Two possible explanations for the observed 
anisotropy were put forward. One is connected to 
the propagation of CR protons expected in direction 
of the nearby spiral arm. The other is that the 
anisotropy might be due to neutron primaries. At 
energies of 10^^ eV neutrons have decay lengths of 
~10 kpc and could thus propagate linearly from 
the GC and Cygnus regions to the position of the 
earth without decaying. More data are needed to 
distinguish between or rule out the models. At 
larger energies the AGASA collaboration did not 
detect any significant large-scale anisotropy with 
respect to the galactic or supergalactic plane [ pO[ . 
The analysed event statistics was 581, 47, and 7 
events with energies larger than 10^^ eV, 4 x 10^^ 
eV, and 10^° eV, respectively. 

The size of the detectors and the time scale 
involved in collecting the data reported on above 
implies that only a new and bold ansatz will allow 
to improve the event statistics in a significant way. 
Especially at the uppermost end of the CR energy 
spectrum where 'New Physics' might be required 
to explain the data, at least an order of magnitude 
improvement in statistics is called for. Currently 
the world statistics of events with energies above 
4 X 10^^ eV is about 200 events corresponding 
to a total integrated exposure of ~1500 km^-yr-sr. 
This energy is called the Greisen-Zatsepin-Kuzmin 
(GZK) cutoff energy and corresponds to the proton 
energy for which the centre-of-mass energy in 
a collision with a cosmic microwave background 
(CMB) photon crosses the pion photoproduction 
threshold, i.e., 



P + 1 



The survival probability of ultrahigh energy CR 
(UHECR) protons of energies 8 x lO^^ eV, lO^o eV, 
and 3 x 10^° eV after having traversed a distance 
of 20 Mpc in the CMB are 0.70, 0.55, and 0.12, 
respectively. Sources of the UHECRs thus have to 
be within the GZK volume in our cosmic vicinity, 
i.e., closer than 30 to 50 Mpc. The possible sources 
discussed in the literature range from nearby proton 
accelerators, i.e., powerful Active Galactic Nuclei 



(AGN) |1[, Gamma-Ray Bursts (GRBs) 
the decay of superheavy relic particles [ p4| , a new 
stable supersymmetric hadron with a mass of a 
few GeV, e.g., the 5*° as a uds-gluino bound state 
which would have a much longer pathlcngth in the 
CMB than ordinary hadrons pq] , to extremely high 
energy neutrinos which would annihilate with relic 
neutrinos to produce hadronic jets |Q . For all these 
types of sources the EHECRs should point back 
(within ~ 5°) to their source(s) if the extragalactic 
magnetic fields are weaker than ^ 0.1 /iG. Should 
these fields be stronger the only sources from above 
for which the directional information is retained 
are the new stable supersymmetric hadron and the 
extremely high energy neutrinos interacting in the 
vicinity of the earth. As recent determinations of 
the magnetic field strength within the local (Virgo) 
supercluster indeed point towards a local field 
strength of ~ 0.5 ^G (e.g., |^) 'conventional' local 
sources like AGN or GRBs might be 'hidden' due 
to large magnetic deflections already over distances 
of the order of Mpc [|8| . 

The experimental investigation of UHECRs will 
be lifted to a new level with the upcoming Pierre 
Auger experiment js^. In March of 1999 the ground 
breaking ceremony for the southern detector took 
place in the Pampa Amarilla in Argentina. The 
detector will consist of ~1600 water Cherenkov 
tanks spread over ~3000 km^ with a 1.5 km grid 
spacing. This air shower detector will be overlooked 
by three fluorescence detectors situated within the 
array. The acceptance of the detector above 5 x 
10^^ eV will be ~14,000 km^-yr-sr. Full operation 
of the detector is planned for the beginning of the 
year 2003. 

4. Gamma Astronomy 

Gamma astronomy is the most recent addition to 
the spectrum of 20th century astronomy. This con- 
cluded the opening of the complete electromagnetic 
spectrum to astronomical observations. Currently 
7-ray astronomy is performed by space-borne and 
ground-based 7-ray detectors in non-overlapping 
energy domains. The space-borne detectors mea- 
sure in the energy domain up to ~10 GeV whereas 
the ground-based detectors are limited to photon 
energies above ^^200 GeV. The remaining gap in 
the electromagnetic spectrum will be closed by the 
upcoming detectors from both sides, i.e., by the 
ground-based and space-borne detectors enlarging 
their energy coverage. The number of 7-ray sources 
detected so far is about 270 in the energy domain 
below 10 GeV |Q and 13 sources (with a varying 

degree of certainty) with E^ ^ 200 GeV il. Note 
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that about 2/3 of the 7-ray sources detected below 
10 GeV have not yet been identified with known as- 
tronomical objects. This class of 'EGRET Uniden- 
tified Sources' thus is one of the main targets for 
upcoming 7-ray experiments. 

Here we only report on some recent results from 
ground-based 7-ray astronomy, i.e., results obtained 
with Imaging Air Cherenkov Telescopes (lACTs). 
The lACT technique is based on the exploitation 
of the inherent differences of 7 and (CR) hadron 
air showers developing in the atmosphere. As a 
function of energy a growing number of secondary 
particles in the air showers have energies in excess 
of the threshold energy of Ethres ~ 21 MeV (at sea 
level) above which Cherenkov radiation is produced 
by the passage of charged particles through air. The 
emission of the Cherenkov light is concentrated into 
a small angular region in the forward direction and 
each shower typically illuminates an area of ^ 5 -10^ 
m^ on the ground. Positioning the I ACT anywhere 
in this area will result in a detection of the incident 
particle. The different images of hadrons and 
photons allow background suppression factors of 
^^600 with stand-alone telescopes. The suppression 
factor is larger by about an order of magnitude 
when one points several (3 to 5) telescopes at 
the same source in the stereoscopic observation 
mode pioneered by the HEGRA collaboration at 
the Canary Island La Palma . Note that due to 
the indirect detection of the 7-rays interferometry, 
however, is not possible. The angular resolution 
of lACTs for individual events is about 0.1 to 
0.2° and the energy resolution is between 10 and 
40%. As for the background suppression, using 
several detectors in stereoscopic mode to measure 
the shower parameters, at current energies (E ~ 
200 GeV) leads to improvements in both parameters 
pl[ . The large effective area gained by the indirect 
detection is one of the major advantages of the 
lACT technique as it allows to detect very low 7-ray 
fluxes. The major drawback of this technique, on 
the other hand, is the low duty cycle. Most lACTs 
are only operated during dark and moonless nights 
resulting in a duty cycle of about 10%. 

Due to the restricted field-of-views of the 
instruments (a few degrees times a few degrees) the 
observations can only be performed in the pointing 
mode. In this mode the telescopes are oriented 
at possible 7-ray sources with guidance from other 
wavelengths together with model predictions for the 
very high energy (VHE) 7-ray channel. One such 
model prediction concerns the important question 
of the origin of the CRs. In most models Supernova 
Remnants (SNRs) are the favoured sites of CR 
acceleration. As the CR nuclei themselves do 



not retain directional information in the galactic 
magnetic field a neutral particle signal must be 
searched for in order to identify the CR sources. 
Whenever hadrons are accelerated to very high 
energies pions should be produced as secondaries. 
However, SNRs are known sources VHE 7-rays 
and the main 7-ray component is due to inverse 
Conipton scattering of ultrarelativistic electrons 
on low energy photons (e.g., CMB photons). In 
order to identify the sources of the CRs one 
thus aims at identifying an additional VHE j-ray 
component stemming from ir^ decay. SNRs thus 
have been intensively studied by the lACTs in the 
past. Although six SNRs have been observed above 
-200 GeV II (Crab nebula, Vela, PSR1706-44, 
SN1006, RJX1713. 7-3946, and (possibly) CasA), 
this additional component has not been clearly 
identified, yet. More sensitive measurements at 
lower energies will be of great importance in 
identifying the spectral component (s) observed in 
the six sources above and to discover 7-ray emission 
in more SNRs. With low energy thresholds, e.g., of 
the upcoming MAGIC Telescope (initially 25 GeV), 
it may be possible to observe a two component j-ray 
spectrum, which should then allow to decouple the 
predicted leptonic and hadronic components in SNR 
shells. In summary, no clear evidence for proton 
acceleration in SNRs was found up to now and the 
sources of the galactic CRs are not identified yet. 

The most spectacular discoveries in ground- 
based 7-ray astronomy are related to the extragalac- 
tic sources of VHE 7-rays. The first such source dis- 
covered in 1992 by the Whipple collaboration [ p3[ 
was the closest BL Lacertae (BL Lac) type object 
Markarian 421 (Mrk 421) at a redshift z = 0.031. 
This AGN at the same time is one of the weakest 
7-ray source detected by the EGRET experiment at 
energies below ~10 GeV Q and showed the short- 
est time variations of its 7-ray emission levels ob- 
served so far with doubling and decaying times in 
the ener gy d omain above —200 GeV as short as 15 
minutes pm . The amplitude ratios between the dis- 
covery fiux level and flux maxima during short out- 
burts have been as large as 30, i.e., between 0.3 and 
10 Crab, where 1 Crab is the (time-independent) 
flux of the standard candle of 7 astronomy, the Crab 
nebula. 

BL Lac objects and fiat spectrum radio-loud 
quasars (FSRQs) are collectively called "blazars". 
The predominantly non-thermal emission shows 
violent variability in most energy bands and 
superluminal motion is observed in VLBI radio 
surveys. The non-thermal emission is believed to 
originate in relativistic jets oriented with small 
angles towards the line of sight. The emission levels 
and time scales observed in the TeV domain, e.g.. 
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for Mrk 421, clearly show that relativistic beaming 
has to be operational in the sources, corroborating 
this model picture. The ultimate energy source is 
believed to be gravitational, i.e., originating from 
accretion of matter onto a supermassive (0(10^ - 
10^) Mq) black hole in the centre of the AGN. How 
the jets are formed or how they are fueled is not 
yet understood. One viable model is based, e.g., 
on the extraction of the rotational energy from the 
ergosphere of the (^maximally) spinning black hole 
p5[ . Irrespective of the actual formation and fueling 
of the jet there are two models put forward to 
explain the observed radiation from 7-ray emitting 
blazars. 

The spectral energy distributions of blazars in 
general seem to consist of two parts. First, a 
low energy component rising from the radio up 
to a broad peak between the infrared (IR) and 
X-ray wavelengths (see, e.g., figure |3|) depending 
on the specific blazar type. This component 
is generally believed to stem from incoherent 
synchrotron emission by relativistic electrons in 
the jets. The origin and composition of the high 
energy component, however, is still a matter of 
debate. Most popular is the model in which the 7- 
rays are produced in optically thin regions through 
inverse Compton scattering of low energy photons 
by the same electron population that produces 
the synchrotron emission. Different models of 
this type differ, e.g., in the origin of the low 
energy seed photons. These photons can be 
synchrotron photons (Synchrotron Self Compton 
Model), photons from the accretion disk, or other 
photons in the vicinity of the jet (External Compton 
Models) (e.g., Q). In a different type of 
model VHE 7-rays are produced by proton-initiated 
cascades (PIC) j9^. In this case normal plasma 
is injected into the jet and shocks within the 
jet accelerate protons and electrons to very high 
energies at which pions can then be photoproduced. 
All or parts of the VHE 7-ray emission in this case 
is the end product of an electromagnetic cascade 
developing in an optically thick acceleration region. 
Note that this model was the only one predicting 
7-ray emission to energies up to and beyond 20 
TeV. A common denominator of all models is that 
the 7-rays originates close to the central black 
hole. Should the low energy fields near to the 
accretion disk, however, be very intense the direct 
escape of the 7-rays becomes difficult due to 77 

e+e" pair production losses VHE 7-ray 

data on time variability and energy spectra taken 
concurrently with data in other energy bands will 
help to constrain or rule out the models. 

The second extragalactic source, Mrk 501, was 
discovered in 1995 [p9|| as a weak source (0.08 Crab). 



In the beginning of 1997 it showed a dramatic 
outburst at TeV energies which lasted for the full 
1997 observation period (March - September) and 
was characterized by wildly varying emission levels. 
The most complete light curve at energies larger 
than 1.5 TeV was obtained by the 6 lACTs of the 
HEGRA collaboration who operated at that time a 
stereoscopic system of 4 lACTs and 2 additional 
stand-alone telescopes. In order to achieve this 
good time coverage of Mrk 501 for the first time 
measurements during moon time were perfo rmed 



with one of the stand-alone telescopes (CTl) [100 



This observation mode considerably increased the 
time coverage with only a moderate increase in 
energy threshold. The light curve measured by 
all HEGRA lACTs is depicted in figure |l| (top 
diagram). Concurrently with HEGRA and many 
other lACTs (e.g., JTol) the AUsky Monitor (ASM) 



of the RXTE X-ray satellite took data on Mrk 501 
in the 2 to 10 keV energy band. The corresponding 
light curve is also shown in figure |l^ (bottom 
diagram). HEGRA performed a correlation analysis 
and found a moderate (0.611 ± 0.057) but clear (8.5 
a) cor relat ion between the TeV and keV emission 
levels |102[ | indicating a common emission region of 
the radiation. 

The determination of the velocity of this 
emitting region within the jet also illustrates the 
impact of such concurrent observations. The bulk 
Lorentz factors, jbuik, of the relativistically moving 
plasma can be determined in two independent 
ways. From radio observations with the Space- 
VLBI technique for Mrk 501 the inclination of 
the jet axis with the line of sight was measured 
to be 6 ~ 10° - 15°. At the same time the 
bulk velocity of the synchrotron emitting plasma, 
f3, was measured to be in the range from 0.990 
to 0.999. These numbers can be translated into 
a Lorentz factor, jtuik, between 7 and 22 or a 
Doppler f acto r S — [ji,uik{^ — P cosQ)]~^ between 
1.3 - 5.6 [ |l03| . The other method is based on the 
simultaneous observation of synchrotron spectra in 
the optical/UV domain and TeV 7-ray spectra and 
makes use of the Synchrotron Self Compton model. 
In this case the concurrent observations require S ~ 
5 1 104 1 consistent with the radio determinations. 

HEGRA investigated the X-ray and TeV 
light curves shown in figure for quasi- 

periodic oscillations (QPOs) ||105|| . In the 
QPO analysis evidence for a 23 day period 
was found in both light curves. In order to 
assign a statistical significance a shot noise model 
assuming as null hypothesis independent, randomly 
distributed flare s in accordance with observations 
was developed [106|. The combined probability 
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Figure 12. Light curves of Mrk 501 in 1997. Top: VHE 7-ray (> 1.5 TeV) emission measured by the HEGRA 
lACTs. Modified Julian Day 50449 correspo nds t o January 1, 1997. Bottom: X-ray light curve measured by the 



RXTE ASM for 2 keV < E < 10 keV. From [102 



for the 23 day period was thus determined to be 
V ~ 2.8 e-04 or 3.5(7. Note that this analysis 
was only possible because the data taken during 
moon time were available. Leaving out these data 
introduces large time gaps and the observed QPO 
signal will be washed out due to the introduction of 
aliasing effects. The interpretation of the observed 
23 day period is complicated by the uncertainty 
of the exact location and environment of the 
emission region. It could be due gravitomagnetic 
precession or g-mode oscillations of the accretion 
disk in analogy to models put forward to explain 
QPOs observed during X-ray flares of the galactic 
analogon to AGN, the X-ray binaries [107|. Or it 
could be due to the action of a close binary system 
of supermassive black holes in the centre of AGN 
which is plausible based on the heavy galaxy merger 
activ ity, e.g., observed in the Hubble Deep Fields 
[108|. Other possibilities are radiating clumps on 
helical trajectories or oblique shock fronts. 

Besides the light curve the other important 
measurement concerns the 7-ray energy spectra 
of AGN during high and low states. From 
their observations of Mrk 501 in 1997 the 
CAT collaboration determ ined energy spectra for 
different 7-ray intensities |109[. Figure O shows 



an example of the multi-wavelength spectra of 
Mrk 501 as measured by the BeppoSAX X-ray 
satellite and the CAT telescope on April 7*'* and 
April 16*'*. The flux in the 7-ray band differs 
by about a factor of 10. In addition CAT found 
that the increase in flux was also accompanied by 
an increase in the spectral hardness |112|. At 
somewhat higher energies this trend could not 
be found in the HEGRA data ||ll||. The most 
striking aspect of the HEGRA energy spectrum of 
Mrk 501 in 1997, besides its independence from 
the flux level, is the maximum energy observed 
from this source. This energy spectrum together 
with the energy spectrum determined for Mrk 421 
from the combined 1997 and 1998 data is shown 
in figure |l|. The 95% CL lower limit on the 
maximum detected energy from this extragalactic 
source is 16 TeV! Due to the distance of Mrk 501 
of about 500 million light years the measurement of 
this maximum energy is a cosmologically important 
signal. The reason for this is that the TeV 
7-rays have to traverse the diffuse extragalactic 
background light (EBL) in order to arrive at the 
location of our Galaxy. The EBL constitutes an 
important cosmological signal as it integrates over 
the emission history of the Universe on Hubble 
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Figure 13. Mrk 501 X-ray and VHE spectra given as 
For April 7"' and 16*, BeppoSAX data were taken 
from lll(| . The EGRET upper limit was taken from 
[111] and corresp ond s to observations between April 9**^ 
and 15*'\ From [|ll2| . 
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Figure 14. The energy spectra of Mrk 421 (1997 and 
1998 combined) and Mrk 501 (199 7) as measured by the 
system of HEGRA lACTs. From fl^]. 



time and length scales. Important parts of the 
EBL spectrum, amongst others the infrared region, 
however, are still only known very poorly from 
direct measurements by satellite detectors. TeV 7- 
rays can interact with infrared photons to produce 
electron-positron pairs. A large density of diffuse 
infrared photons therefore results in measurable 
absorption features in extragalactic 7-ray spectra. 
As the energy dependence of the EBL flux is roughly 
oc quasi- exponential cutoffs of the TeV spectra 
are expected. By turning this argument around 
one sees that a beam of TeV 7-rays traveling 
cosmological distances can be used to probe this 
background field. As the relevant Thomson cross 
section is strongly peaked at the threshold this 
probing can even be performed spectroscopically by 
observing sources at different distances. Because 
of the energy dependence of the EBL spectrum 
the cutoff condition, i.e., optical depth reaching 
unity, for sources at different cosmological distances 
occurs at different 7-ray energies. A 'cosmological 
7-ray horizon with the Universe opening up as one 
goes towards lower 7-ray energies should result. A 
small fraction of the EBL spectrum is shown in 
figure [ij. What is striking is that in the shown 
energy interval between the far infrared ('^200 fim) 
and the ultraviolet (UV) (~0.1 /im) the figure is 
dominated by wildly scattered upper and lower 
limits. The main reason for this is the enormous 
foreground radiation in this spectral range by the 
zodiacal light in our solar system. Besides the 
limits and measurements obtained with 'direct' 
methods also the upper limits on the IR EBL level 
derived from the TeV 7-ray spectrum of Mrk 501 
as measured by the Whipple collaboration is shown 
in figure |l^. A similar analysis was performed on 
the basis of the HEGRA data resulting in an upper 
limit of 1.1 X 10^'^ eV/cm^ at ^25 /im assuming a 
specific shape of the IR background . As visilale 
in figure ^the application of the methods of TeV 7- 
ray astronomy have already led to an improvement 
of the upper limits in the mid-infrared range by 
more than an order of magnitude. 

Especially in the mid-infrared range where 
the models predict a pronounced dip in the 
diffuse extragalactic photon flux the dominating 
foreground will also in the future probably not 
allow a certain direct determination of the flux 
level by infrared satellites. The level of the EBL 
in this range, however, is an important signal of 
structure formation in the Universe and constitutes 
a (complicated) convolution of star formation rate, 
initial mass function and its evolution, and the dust 
history of the Universe. It is therefore of paramount 
interest to determine the flux with the methods of 
TeV astronomy. As this means measuring cutoff 
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Figure 15. The density of the EBL between the far- 
infrared and the ultraviolet. Filled squares with arrows 
indicate 95%CL upper limits by various experiments. 
Open squares are detections (140 and 240 /^m) and a 
tentative detection (3.6 ^m, denoted by T) by DIRBE. 
The filled circles are lower limits from IRAS and ISO 
galaxy counts. For references see [114|. The horizontal 
bars with arrows indicate the u pper limits from the Mrk 
501 spectrum derived in |115|. The dashed line above 
300 ^im indicates the level of the 2.7K CMB. The dotted 
and dot-dashed hues are model predictions from jll(j ] 
and the dashed line a prediction taken from [117|. From 



energies in the TeV domain, it makes it necessary 
to find more close-by (on cosmological scales) TeV 
sources, i.e., raising the instruments sensitivities as 
planned for the next generation of 7-ray detectors. 

The final result from TeV 7-ray astronomy 
reported on here concerns the test of Lorentz 
invariance made possible by the combination of high 
energy emission and cosmological distances. Some 
of th e ansatze for Quantum Gravity (QG), e.g., 
1 118 1, yield non-Lorcntz- invariant terms which lead 
to modified laws of propagation and interaction of 
neutral particles as a result of interactions with 
the quantum gravity medium. Measurable time 
delays can be expected if the particles have energies 
close to the QG scale (expected to be close to 
the Planck scale, i.e., 10^^ GeV) or if they have 
traversed cosmological distances. Biller et al. |119 



analysed the shortest flare of Mrk 421 observed on 
May 15 1996 by the Whipple collaboration Q for 
a time delay between 7-rays with energies less than 
1 TeV and those with energies above 2 TeV. No 
difference was found within the measurable shortest 
time bin of 280 s and a lower limit of 4 x 10"'^^ 
GeV could be placed on the relevant QG energy 



scale. In general this type of investigation can be 
used to place limits on any non-Lorcntz invariant 
propagation term. By searching for effects of this 
kind in very intense transient phenomena occurring 
at cosmological distances like, e.g., GRBs, the 
upcoming MAGI C te lescope will have sensitivity at 
the Planck scale |12C|. For comparison, the current 
best limits are about 1% of the Planck scale. 

5. Neutrino Astronomy 

The observation of low energy astrophysical 
neutrinos from the sun and SN1987A has led 
to important astrophysical and particle physics 
results. The fundamentally new aspect of 
performing neutrino astronomy in contrast to 
astronomy with electromagnetic radiation lies in the 
penetration power of neutrinos, i.e., also optically 
thick regions can be studied in the neutrino 'light'. 
In addition neutrinos do not suffer attenuation 
losses due to absorption on background radiations 
as do 7-rays or ultrahigh energy hadrons. We know 
from the observation of CRs that large numbers 
of hadrons are accelerated to very high energies 
in astrophysical environments. Whenever processes 
accelerate hadrons to VHE energies besides the 
production of 7-rays we naturally expect neutrino 
ffuxes as least as large as the 7-ray ffuxes and 
neutrino and 7-ray astronomy are complementary 
approaches to many astroparticle physics questions. 
Also several astrophysical events, e.g., supernovae 
(SN) or GRBs, are predicted to emit their maximum 
power in the neutrino channel. In summary, all 
environments where photoproduction of pions is 
likely to occur are potential high energy neutrino 
sources. These are the (unknown) sources of 
CRs, the jets of AGN if protons are accelerated 
within the jets, the Galactic disk because of CR 
interactions with the interstellar matter, or the 
centres of galaxy clusters where UHECRs might 
be present. In addition there are exotic sources 
like the pairwise annihilation of neutralinos and 
radiation from topological defects. Note that the 
detection of AGN as point sources for high energy 
neutrinos is considered to be the experimentum 
crucis in distinguishing between jet models where 
only electrons and positrons or where 'normal 
plasma', i.e., electrons and protons produce the 
observed VHE 7 radiation. Although the neutrino- 
proton cross section rises linearly up to energies of 
~ 10^ GeV before ffattening due to W propagator 
effects the strongly falling fluxes expected from 
most conceivable sources demands that high energy 
neutrino detectors have to be huge instruments. For 
more information a nd the p otential of this exciting 
new field see, e.g., [122, 123(|. 
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The instruments that are taking data 
(BAIKAL, AMANDA-B13), are under construction 
(AMANDA-B, ANTARES Phase II), or are in the 
test phase (NESTOR) are huge detectors with 
effective collection areas between 2 xlO^ and ~10^ 
m^. Due to the expected low rates, however, all 
these detectors are considered to be first or second 
steps towards the construction of one or two 1 km'^ 
detectors (e.g., ANTARES, ICECUBE, NESTOR). 
Basically the detectors are muon detectors aimed 
at measuring i/^ and induced upward going 
high energy muons. At high energies this gives the 
added advantage of enlarging the fiducial volume of 
the detectors due to the long range of high energy 
muons in water or ice. Water or ice are the detector 
media of choice because the detection of the muons 
utilizes the Cherenkov light emitted by the water 
or ice when traversed by relativistic particles with 
energies above the Cherenkov threshold. The main 
parameters of these neutrino telescopes are the 
angular and energy resolution. In general water 
detectors will have a better angular resolution but 
worse energy resolution compared to ice detectors. 
In addition and Uf, fluxes can be easily measured 
at large energies based on the detection of induced 
electromagnetic showers. The major background to 
extraterrestrial neutrinos is given by the downward 
going atmospheric muons and the upward going 
muon flux induced by atmospheric neutrinos both 
of which are produced by the interaction of CRs in 
the earths atmosphere. In addition to constituting 
two sources of background, however, these two 
beams can be used as test beams for proving the 
feasibility of these detectors by measuring the 
fluxes and angular distributions of these beams at 
very high energies. 

An ice detector currently under construction is 
the AMANDA-B detector at the South Pole. It 
will have an effective collection area of 10,000 m^, 
an energy threshold of ~50 GeV and an angular 
resolution of about 2.5° per muon track [123|. While 
being constructed at different depths the AMANDA 
detector has been taking data. The AMANDA- 
A detector was installed in a depth between 810 
and 1000 m below the south pole's surface and is 
characterized by a high concentration of residual air 
bubbles leading to strong scattering of light. While 
this dramatically worsened the angular resolution 
it improved the calorimetric performance and the 
1995 data cascade trigger data of AMANDA-A 
could be used to put an upper limit on the diffuse 
flux o f Ve. + hi the energy range from 10^ to 10^ 
GeV (l2|l. Due to the small size of AMANDA- 
A this limit is, however, not yet restricting the 
model space of fluxes predicted in AGN models 
(see figure |l6|). After AMANDA-A the installation 
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Figure 16. Limits on the diff use Ve -\- Ve flux 
as determined b y th e Frejus |125l and AMANDA-A 
detectors. From [124|. 



towards the AMANDA-B detector started with the 
deployment of 4 strings at a depth between 1.5 
and 2 km. For results from this detector see 
126|. The data taken during the first year of 



operation of the 10 string AMANDA-BIO detector 
installed in the depth between 1500 and 1980 m 
were analyzed for upward going muon candidates 
[ |127| . In a preliminary analysis of data taken during 
113 days 17 events passed all selection cuts. This 
corresponds to a reduction of the number of events 
in the solid angle where the 'pencil like' AMANDA- 
BIO detector is sensitive by more than a factor of 
10^. The result was found to be in accordance with 
the Monte Carlo expectation of 21.1 events. The 
reconstructed zenith angle distributions for the full 
data set (4.9x10* events), and the data sets after 
different selection cuts are shown in figure |l^. 

The BAIKAL detector, operational since April 
1993, also took data while being enlarged to its 
final NT-200 configuration in April 1998. BAIKAL 
can be considered to be the water detector which 
proved the feasibility of very high energy neutrino 
telescopes when it reconstructed the first two 
atmospheric neutrino candidates from the initial 
NT-36 data taken between April 1994 and March 
1995 [^2|]. In its NT-96 configuration it took 
data from April 1996 until March 1997 and a very 
tight limit could be placed on the diffuse 



fiux for neutrino energies larger than 10 TeV |129|. 
This limit together with some other experimental 
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Figure 17. Reconstructed zenith angle distributions 
of 113 days of AMANDA-BIO data after quality cuts 
tightening from top to bottom. The 17 events with 
a cos9 zenith clos e to -1 are atmospheric neutrino 



candidates. From |127 



limits and model predictions is shown in figure 
Although the BAIKAL NT-96 limit is still about 
a factor of 10 above the highest model predictions 
it demonstrates the capability of the underwater 
technique. This technique is also going to be applied 
by the ANTARES project planning the construction 
of a 1 km'^ detector in the Mediterranean Sea (e.g., 
[13l[). Since 1997 ANTARES is conducting an 
intense R&D program aimed at the installation of 
a 0.1 km^ detector in ANTARES Phase II by 2002 
to 2003. In this approved phase 13 strings with 
- 1000 PMTs will be deployed in about 2500 m 
depth. Monte Carlo simulations of the performance 
look very promising, e.g., an angular resolution of 
~0.2° is indicated. 

6. Outlook 

The major potential of Astroparticle Physics is 
given by the synergistic effects of combining par- 
ticle, nuclear and atomic physics with astrophysics 
and cosmology and it has started to make an impact 
on important physics questions. The improvement 
of detector technology coupled with improved 'cali- 
brations' of astrophysical beams, e.g., CR fluxes or 
7 and v fluxes from extragalactic objects, will lead 
to an even larger sensitivity to fundamental ques- 
tions in the future. 
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